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We present molecular dynamics simulations of the motion of a single rigid rod in a disordered 
static 2d-array of disk-like obstacles. Two different configurations have been used for the latter: 
A completely random one, and which thus has an inhomogeneous structure, and an homogeneous 
"glassy" one, obtained from freezing a liquid of soft disks in equilibrium. Small differences are 
observed between both structures for the translational dynamics of the rod center-of-mass. In 
contrast to this, the rotational dynamics in the glassy host medium is strongly slowed down in 
comparison with the random one. We calculate angular correlation functions for a wide range 
of rod length L and density of obstacles p as control parameters. A two-step decay is observed 
for large values of L and p, in analogy with supercooled liquids at temperature close to the glass 
transition. In agreement with the prediction of the Mode Coupling Theory, a time- length and 
time-density scaling is obtained. In order to get insight on the relation between the heterogeneity of 
the dynamics and the structure of the host medium, we determine the deviations from Gaussianity 
at different length scales. Strong deviations are obtained even at spatial scales much larger than the 
rod length. The magnitude of these deviations is independent of the nature of the host medium. 
This result suggests that the large scale translational dynamics of the rod is affected only weakly 
by the presence of inhomogeneities in the host medium. 

Ref.: AIP Conference Proceedings 708 (2004) 576-582 



I. INTRODUCTION 

Since it was initially introduced by Lorentz as a model 
for the electrical conductivity in metals 0|, the prob- 
lem of the Lorentz gas has given rise to a substan- 
tial theoretical effort aimed to understand its proper- 
ties H,H,^|MIEI^H]' In this model, a single classical 
particle moves through a disordered array of static ob- 
stacles. It can thus be used as a simplified picture of the 
motion of a light atom in a disordered environment of 
heavy particles having a much slower dynamics. In the 
simplest case, where the diffusing particle and the obsta- 
cles are modeled as hard spheres, an exact solution for 
the diffusion constant exists in the limit of low densities 
of obstacles 0. However, the problem becomes highly 
non-trivial with increasing density, where dynamic cor- 
relations and memory effects start to become important 
for the motion of the diffusing particle, and the system 
shows the typical features of the dynamics of supercooled 
liquids or dense colloidal systems, such as a transition to 
a non-ergodic phase of zero diffusivity In par- 

ticular, diffusion constants and correlation functions are 
non-analytical functions of the density 0) S S • More- 
over, correlation functions show non-exponential long- 
time decays. 

Diffusing particles and obstacles are generally modeled 
as disks or spheres in two and three dimensions, respec- 
tively. Much less attention has been paid to systems that 
have orientational degrees of freedom. Motivated by this 
latter question, we have recently started an investigation, 
at low and moderate densities of obstacles, on a general- 
ization of the Lorentz gas, namely a model in which the 
diffusing particle is a rigid rod • An array of randomly 



distributed disks has been used for the host medium. For 
simplicity simulations have been done in two dimensions, 
reducing the degrees of freedom to the center-of-mass po- 
sition and the orientation of the rod axis. As in the case 
of supercooled liquids llll 113 . Il5| or dense col- 

loidal systems 0, 0, Il8l| . one observes at intermediate 
times a caging regime for the rod center-of-mass motion, 
which is due to the steric hindrance produced by the 
presence of the neighboring obstacles. 

More interestingly, strong deviations from Gaussian- 
ity have been obtained for the incoherent intermediate 
scattering function at wavelengths much longer than the 
rod length, giving evidence for a strongly heterogeneous 
character of the long-time dynamics at such length scales. 
The inhomogeneous structure of the model used for the 
host medium has been pointed out as a possible origin of 
such non-Gaussian effects since a random configuration 
features large holes on one side but on the other side also 
dense clusters of obstacles. The presence of holes might 
lead to a finite probability of jumps that are much longer 
than the average "jump length" and thus to a heteroge- 
neous dynamics. 

In order to shed new light on this question, we have re- 
peated the simulations by taking a disordered but homo- 
geneous configuration of the obstacles instead of a ran- 
dom one. We will see, however, that large scale non- 
Gaussian effects are not significantly affected by the par- 
ticular choice of the configuration of the obstacles. Some 
information on angular correlation functions is also pre- 
sented. The article is organized as follows: In Section II 
we present the model and give details of the simulation. 
Translational and angular mean-squared displacements 
are presented in Section III. The behavior of angular cor- 
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relation functions is shown and discussed in Section IV, 
in terms of the Mode Couphng Theory. Non-Gaussian ef- 
fects in the random and glassy models of the host medium 
are investigated in Section V. Conclusions are given in 
Section VI. 



II. MODEL AND DETAILS OF THE 
SIMULATION 

The rigid rod, of mass M, was modeled as N aligned 
point particles of equal mass m = M/N, with a bond 
length 2a. The rod length is therefore given by L = 
{2N — l)cr. The positions of the obstacles in the ran- 
dom configuration were generated by a standard Pois- 
son process. In order to obtain the glassy host medium, 
we equilibrated at a reduced particle density p — 0.77 
and at temperature T — e/fce a two-dimensional ar- 
ray of point particles interacting via a soft-disk potential 
V{r) = e((T/r)^^. This procedure produced an homoge- 
neous liquid-like configuration. The latter was then per- 
manently frozen and was expanded or shrunk to obtain 
the desired density of obstacles, defined as p = «obs/^box' 
with riobs the number of obstacles and Zbox the length 
of the square simulation box used for periodic boundary 
conditions. 

The same soft-disk potential V{r) was used for the 
interaction between the particles forming the rod and 
the obstacles. For computational efficiency, V{r) was 
truncated and shifted at a cutoff distance of 2.5(7. In 
the following, space and time will be measured in the 
reduced units a and (cr^m/e)^/^, respectively. Typically 
600-1000 realizations of the ensemble rod-obstacles were 
considered. The set of rods was equilibrated at T = 
e//cB- After the equilibration, a production run was done 
at constant energy and results were averaged over the 
different realizations. These runs covered 10^ time units, 
corresponding to (1 — 5) • 10^ time steps, depending on the 
step size used for the different rod lengths and densities. 
Run times were significantly longer than the relaxation 
times of the system. 



III. MEAN-SQUARED DISPLACEMENTS 

Figs, la and lb show respectively for p = 6 • 10^'^ 
a comparison between the random and the glassy con- 
figuration for the mean-squared displacement of the rod 
center-of-mass ((Ar(t))2) and the mean-squared angular 
displacement ((A$(i))^). Brackets denote ensemble av- 
erage. In order to facilitate the observation of the differ- 
ent dynamic regimes, data have been divided by the time 
t. The comparison covers a wide range of rod lengths 
from L ~ O.ldnn to L ~ lOdnn, where d^n = p^^^^ ~ 13, 
is the average distance between obstacles for the men- 
tioned density. At short times the rod does not feel the 
presence of the neighboring obstacles and ((Ar(t))^) and 
((A$(<))^) show the quadratic time-dependence charac- 
teristic of a ballistic motion. For the shortest rods a sharp 
transition to the long-time linear regime is observed. In 
contrast to this, for long rods a crossover regime between 
both limits, showing a weaker time dependence that the 




10"^ 10"' lo" lo' 10^ 10^ lO"* 10^ t 10^ 



FIG. 1: Mean-squared displacement of the center-of-mass (a) 
and mean-squared angular displacement (b), both divided by 
the time t, for p = 6 • 10""^ and different rod lengths, for the 
random (solid lines) and glassy (dashed lines) configuration 
of the obstacles. 



ballistic motion, is present over 1-2 decades of interme- 
diate times. Such a crossover corres pon ds to the well- 
known caging regime |il[lllil[ll(ll[lllll[li ob- 
served in supercooled liquids or dense colloidal systems. 
Due to the presence of the neighboring obstacles, the par- 
ticle is trapped within an "effective cage" for some time 
until it escapes from it and begins to show a diffusive 
behavior. 

No significant differences are observed between the val- 
ues of ((Ar(t))^) for a same rod length in the two dif- 
ferent, random and glassy, configurations. This is more 
clearly seen by calculating for each L the ratio IJJ^m /-^cm 
between the center-of-mass translational diffusion con- 
stant 13cm in both configurations. The latter is calcu- 
lated as the long-time limit of ((Ar(i))^)/4t. As shown 
in Fig. 2, differences in Dqm between the random and 
the glassy host medium are less than a factor 1.5. The 
values of Dqu are systematically lower for the glassy 
host medium, evidencing a stronger backscattering for 
the diffusion, as intuitively expected from the homoge- 
neous character of this latter configuration in contrast to 
the random one (see Fig. 3), where the presence of holes 
facilitates diffusion. It is noteworthy that the maximum 
difference between the values of Dcm in both configu- 
rations takes place for L ^ 10 — 20, i.e., when the rod 
becomes longer than dnn- Thus, while in the homoge- 
neous glassy host medium the rod will not be able to 
pass transversally between two neighboring obstacles, in 
the inhomogeneous random configuration this diffusion 
channel will still be present due to the presence of holes 
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FIG. 2: Ratio between the rotational and center-of-mass 
translational difltusion constants for the random and the 
glassy host medium at p = 6 • 10~^ and different rod lengths. 




FIG. 3: Two configurations of the obstacles for a same den- 
sity p = 6 • 10""^. Left: A random configuration. Right: 
A homogeneous glassy configuration. Short and long arrows 
correspond, respectively, to distances of 40 and 100. 

and will only be suppressed for very long rods. 

Strong differences are observed between both mod- 
els of the host medium in the case of the rotational 
dynamics, as shown in Fig. lb for ((A$(i))^} and in 
Fig. 2 for the ratio of the rotational diffusion constants, 
Dj^'^/Df^'^. is calculated as the long-time limit of 

((A$(t))^}/2i.) For rod lengths smaller than dnn the ro- 
tational dynamics is only weakly sensitive to the config- 
uration of the medium and the ratio D'^'^/D^^ remains 
close to unity. However, for rods longer than dnn the lat- 
ter ratio strongly increases up to a maximum of ~ 4.5 for 
L ^ 40. The position of the maximum can be again ratio- 
nalized from the inhomogeneous structure of the random 
host medium. Short arrows in the configurations of Fig. 
3 for p = 6 • 10^'^ indicate distances of 40. From the com- 
parison between both configurations it is clear that while 
in the glassy medium rods of L 40 can perform only 
small rotations within the tubes formed by the neighbor- 
ing obstacles, in the random medium they can go into the 
holes and thus rotate freely over a larger angle. There- 
fore, the long-time angular displacement will grow up 
more quickly than for the motion between narrow tubes 
in the glassy configuration. It must also be mentioned 
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FIG. 4: Correlation function (cos mA<l>(t)) for L = 99 and 
p — 6 - W~^, with m = 1, 2, 3. ..12 (from top to bottom). Data 
correspond to the random host medium. 

that even for the largest investigated rod length, L — 99, 
where motion between tubes also dominates the diffusion 
in the random medium, the ratio D^'^/D^^ is still sig- 
nificantly different from unity. Thus, while in the glassy 
medium the walls of the tube are formed by uniformly 
distributed obstacles, in the random configuration long 
distances are allowed between some of the neighboring 
obstacles forming the tube, leading to additional escap- 
ing channels (see long arrows in Fig. 3). In principle, 
only for extremely long rods the ratio Dj^'^/D^'^ is ex- 
pected to approach unity. 



IV. CORRELATION FUNCTIONS 

The Mode Coupling Theory (MCT) [H IH IH HI 
makes precise quantitative predictions on the dynamics 
of supercooled liquids or dense colloidal systems. In its 
idealized version, it predicts a dynamic transition from 
an ergodic to a non-ergodic phase at some critical value 
of the control parameters. These are usually the tem- 
perature or the density, though in principle the MCT 
formalism can be generalized to other control parame- 
ters. Moreover, it has been recently tested in a kind of 
systems very different from liquids or colloids, such as 
plastic crystals suggesting a more universal charac- 
ter for this theory. Motivated by this possibilities, we test 
some predictions of MCT for the rotational dynamics of 
the rod with L and p as control parameters. 

According to MCT, upon approaching the critical 
point, the correlation function shows a two-step decay. 
The initial decay corresponds to the dynamic transition 
from the ballistic to the caging regime. As a consequence 
of the temporary trapping in the cage formed by the 
neighboring obstacles, the particle does not lose the mem- 
ory of its initial position and the correlation functions de- 
cay very slowly, giving rise to a plateau at intermediate 
times between the ballistic and the diffusive regime. The 
dynamics close to this plateau is usually referred as the 
/3-relaxation. The closer the control parameters are to 
the critical values, the slower is the mobility of the par- 
ticles and the longer is the lifetime of the cage, leading 
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FIG. 5: (a) Correlation function (cos3A$(t)) for p = 6-10"^ 
and different rod lengths; (b) The same function for L = 99 
and different densities. Solid and dashed lines correspond, 
respectively, to the random and the glassy configurations of 
the host medium. 



FIG. 6: (a) Time-length superposition for (cos A<l?(f)) at den- 
sity p — 6 ■ 10~^ ; (b) Time-density scaling for the same func- 
tion for length L = 99. For the sake of clarity the lowest 
values of p and L have been represented with points. Dashed 
lines are fits to KWW functions. The stretching exponents 
are indicated in the figure. Data correspond to the random 
host medium. 



to a longer plateau in the correlation function. Finally, 
at much longer times, the particle escapes from the cage 
and loses memory of its initial environment, leading to 
a second long-time decay of the correlation function to 
zero, known as the a-relaxation. 

Fig. 4 shows the behavior of the angular correlators 
(cosmA$(i)) for m = 1,2,3, ...12 for a density of obsta- 
cles p = 6 • 10^"^ and a rod length L — 99 much longer 
than dnn. Data are shown for the random host medium, 
though for the glassy one they show the same qualitative 
behavior. Though it is difficult to see it for the small- 
est values of to, the plateau is clearly visible for m > 4. 
Thus, long rods will only perform small rotations be- 
fore hitting the walls of the tube. As a consequence, the 
correlators for small to will decrease only very weakly 
during the ballistic regime, and hence it will be diffi- 
cult to see the subsequent plateau. In contrast to this, 
for higher-order correlators, the angle rotated during the 
ballistic regime will be amplified by a factor to, leading 
to a stronger decay before the caging regime and facili- 
tating the observation of the plateau. The latter begins 
to develop around t ~ 30. As can be seen in Fig. lb for 
the curve for L — 99, this time corresponds to the begin- 
ning of the caging regime for the mean-squared angular 
displacement, in agreement with the MCT prediction. 

Fig. 5 shows for to = 3, and for the two models of 
host medium, how the plateau develops with increasing 
rod length and density of obstacles. While for values of L 
below dnm the correlators show a simple decay and within 
the error bar show no difference between both models of 
the host medium, a clear difference is observed for longer 



rods. Thus, for the glassy medium rotational relaxation 
times become about a decade larger and the plateau is 
significantly higher than in the random configuration, in 
agreement with the results presented in Section III for 
the mean-squared angular displacement. 

Another important prediction of the MCT is the so 
called "second universality": given a correlator C)j 
with ^ a control parameter, then one has that, in the 
time scale of the a-regime, the correlator follows a scaling 
law g{t,C,) = g{t/T{C,)), where t{() is the C-dependence 
of the relaxation time r of the a-regime for such a corre- 
lator, and g is a master function. The relaxation time is 
in practice defined as the time where the correlator de- 
cays to an arbitrary but small fraction of its initial value, 
or is obtained from fitting the a-decay of the correla- 
tor to a Kohlrausch- Williams- Watts (KWW) function, 
exp( — (i/r)'^), widely used in the analysis of relaxations 
in complex systems l24j . In most of the experimental sit- 
uations the relevant control parameter is the temperature 
and for that reason, the second universality is often re- 
ferred as the "time-temperature superposition principle" . 
Now we test the existence of a time-length and a time- 
density superposition principle for the angular correlators 
of the rod. Fig. 6 shows the correlators (cos A$(t)) as a 
function of the scaled times t/T{L) for constant density 
p — 6 ■ 10~^, and t/T{p) for constant rod length L — 99. 
Data are shown for the random host medium. The re- 
laxation times T have been defined as (cos A<I>(t)) ~ 1/e. 
Apart from the trivial scaling in the limit of short rods 
and low densities, a good superposition to a master curve 
is obtained for larger values of L and p, confirming the 



second universality of the MCT for this system. The 
master curve can be well reproduced by a KWW func- 
tion with a stretching exponent [3 — 0.8, as shown in 
Fig. 6. Analogous results, with a very similar stretching 
exponent, are obtained for the glassy host medium. 



V. NON-GAUSSIAN EFFECTS 

In Einstein's random walk model for diffusion, par- 
ticles move under the effect of collisions with the oth- 
ers. When a particle undergoes a collision, it changes 
its direction randomly, and completely loses the mem- 
ory of its previous history. When the time and spatial 
observational scales are much larger than the charac- 
teristic time and the mean free path between collisions, 
the van Hove self-correlation function, i.e., the time and 
spatial probability distribution Gs(r,t) of a particle ini- 
tially located at the origin, is given by a Gaussian func- 
tion Jisl. This functional form is exact for an ideal 
gas and for an harmonic crystal. It is also valid in the 
limit of short times, where atoms behave as free par- 
ticles. When the system is observed at time and length 
scales comparable to those characteristic of the collisions, 
the possible intrinsic dynamic heterogeneities of the sys- 
tem will be reflected by strong deviations from Gaus- 
sianity in Gg (r, t) . Such deviations are usually quan- 
tified by the so-called second-order non-Gaussian pa- 
rameter Q!2(i), which in two-dimensions is defined as 
a2{t) = [((Ar(t))4)/2((Ar(t))2)2] _ l. For a Gaussian 
function in two dimensions, a2{t) = 0, while deviations 
from Gaussianity result in finite values of a2{t). 

Our previous investigation [1(TI | on the non-Gaussian 
parameter in the random host medium revealed some fea- 
tures similar to those observed in supercooled liquids or 
dense colloids, such as the development of a peak (see 
also Fig. 7), which grows with increasing rod length - 
in analogy to decreasing temperature in supercooled liq- 
uids or increasing density in colloids. As also observed in 
these latter systems [H iH El ElIIl III 113 El , the re- 
gion around the maximum of the peak corresponds to the 
time interval corresponding to the end of the caging and 
the beginning of the long-time diffusive regime. Thus, 
the breaking of the cage leads to a finite probability of 
jumps much longer that the size of the cage, resulting in 
a strongly heterogeneous dynamics at that intermediate 
time and spatial scale, which is reflected by a peak in the 
non-Gaussian parameter. 

Another interesting result was the observation that, in 
particular for long rods, in the time span of the simula- 
tion, i.e. for time scales that are much longer than the 
typical relaxation time of the system, a2{t) did not de- 
cay to zero but remained finite [lOj. Thus, the long-time 
dynamics is still significantly non-Gaussian, i.e. hetero- 
geneous, at the spatial scale -much larger than the rod 
length- covered by the simulation. In order to investi- 
gate the possibility of a relation of this large-scale non- 
Gaussian dynamics with the presence of holes in the ran- 
dom structure of obstacles, that might lead to a finite 
probability of jumps much longer than the average, as 
pointed out in Ref.[10], we have calculated Q2(i) for the 
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FIG. 7: Time-dependence of the non-Gaussian parameter 
0^2 (t) for different values of L and for the random and the 
glassy host medium at density p — Q ■ 10~^. 

dynamics of the rod in the glassy host medium. Fig. 7 
shows a comparison between the non-Gaussian parame- 
ter in both structures for several rod lengths and den- 
sity p = & ■ 10""^. While in the limit of short and long 
rods a2{t) takes similar values for both configurations 
in all the time window investigated, a much more pro- 
nounced peak is observed at intermediate rod lengths for 
the glassy configuration. This result can be rationalized 
by the absence of holes in the homogeneous glassy host 
medium. Thus in the latter, jumps that lead to an escape 
from the cage formed by the neighboring obstacles will 
be necessarily long, since long rods will be confined in 
tubes which they will be able to leave only by making a 
long longitudinal motion. In contrast to this, the above 
mentioned inhomogeneous nature of the tube walls in the 
random structure will allow the rod to escape the cage 
also by shorter jumps, resulting in a less heterogenous 
caging regime. Only rods much longer than the hole size 
will need long jumps to escape from the tubes in the ran- 
dom host medium and the caging regime will become as 
heterogeneous as in the glassy host medium, as evidenced 
by the similar peak heights in Fig. 7 for L = 59 and 99. 

Concerning the long-time dynamics, the non-Gaussian 
parameter also remains finite for the glassy host medium, 
and no important differences with the random configura- 
tion are observed. Therefore, contrary to what was pre- 
viously pointed out El|i the non-Gaussianity observed at 
large length scales for the long-time dynamics in the in- 
homogeneous random host medium is also present in the 
homogeneous glassy one, and is not related to the pres- 
ence of holes in the configuration of obstacles, the latter 
having effects only in the intermediate caging regime. 

Another way of visualize the non-Gaussian effects 
at large length scales is obtained by representing the 
intermediate incoherent scattering function Fs{q,t) = 
(exp(— jq • Ar(t))), at very long wavelengths A ~ '2,'K/q. 
Brackets denote ensemble and angular average. Fig. 8 
shows this latter function in the random host medium 
for L = 99 and p — & ■ 10~^, at different (long) wave- 
lengths. (Note that the simulations have been extended 
one order of magnitude respect to those presented in 
Ref. [Hi.) The presence of non-Gaussian effects are made 
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FIG. 8: Intermediate incoherent scattering function Fs{q = 
2Tv/\,t) at different wavelengths A for p = 6 ■ 10~^ and 
L = 99. Points correspond to simulation data in the ran- 
dom host medium. Lines are the curves in Gaussian approx- 
imation as calculated from the center-of-mass mean-squared 
displacement obtained from the simulation (see text). 

clear by comparing the different curves with those cor- 
responding to the Gaussian case in two-dimensions [isf . 
exp[— ((Ar(t))^)q^/4], where we use the center-of-mass 
mean-squared displacement ((Ar(t))^) calculated from 
the simulations. From this comparison it is clear that 
significant non-Gaussian effects are still present at length 
scales of at least ~ 9000, i.e. two orders of magnitude 
larger than the rod length. Whether this result indicates 
that the diffusion of a rod in a disordered array of obsta- 
cles is actually a Poisson process at any length scale is 
an open question. 

VI. CONCLUSIONS 

By means of molecular dynamics simulations, we have 
compared the dynamics of a rigid rod in two models of a 
2d-disordered static host medium: a random configura- 
tion of soft disks and another "glassy" one obtained from 
freezing a liquid of soft disks in equilibrium. While the 
former is characterized by the presence of big holes and 
clusters of close obstacles, the latter presents an homo- 
geneous structure. No significant differences have been 
observed in the translational dynamics of the rod center- 
of-mass. However, rotations are much more hindered in 
the glassy host medium. 

Angular correlation functions have been calculated for 
a wide range of rod length and density of obstacles. In 
agreement with the predictions of the Mode Coupling 
Theory, these functions show a plateau at the time scale 
of the caging regime, and follow a time-length and a time- 



density scaling for the long-time dynamics. 

Strong non-Gaussian behavior has been observed at 
large length scales, though no significant differences are 
evidenced between the random and the glassy configu- 
ration of the obstacles. This result suggests that the 
long-time translational dynamics is not controlled by the 
presence of inhomogeneities in the host medium. 
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ropean Community's Human Potential Program under 
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